Saturation of azimuthal anisotropy in Au + Au collisions at y/s^^ = 62 - 200 GeV 



S.S. Adler, 5 S. Afanasiev, 20 C. Aidala, 5 - 10 N.N. Ajitanand, 48 Y. Akiba, 23 ' 42 ' 43 A. Al-Jamel, 37 J. Alexander, 48 
R. Amirikas, 14 K. Aoki, 27 ' 42 L. Aphecetche, 50 R. Armendariz, 37 S.H. Aronson, 5 R. Averbeck, 49 T.C. Awes, 38 
B. Azmoun, 5 R. Azmoun 49 V. Babintsev, 17 A. Baldisseri, 11 K.N. Barish, 6 P.D. Barnes, 30 B. Bassalleck, 36 
S. Bathe, 6 ' 33 S. Batsouli, 10 V. Baublis, 41 F. Bauer, 6 A. Bazilcvsky, 5 - 17 ' 43 S. Belikov, 5 ' 17 - 19 R. Bennett, 49 
Y. Berdnikov 45 S. Bhagavatula, 19 M.T. Bjorndal, 10 J.G. Boissevain, 30 H. Borel, 11 S. Borenstein, 28 K. Boyle 49 
M.L. Brooks, 30 D.S. Brown, 37 N. Bruner, 36 D. Bucher, 33 H. Buesching, 5 ' 33 V. Bumazhnov, 17 G. Bunce, 5 - 43 
J.M. Burward-Hoy, 29 * 30 ' 49 S. Butsyk, 49 X. Camard, 50 S. Campbell, 49 J.-S. Chai, 21 P. Chand, 4 W.C. Chang, 2 
S. Chcrnichcnko, 17 CY. Chi, 10 J. Chiba, 23 M. Chiu, 10 I.J. Choi, 57 J. Choi, 22 R.K. Choudhury, 4 T. Chujo, 5 ' 54 
V. Cianciolo, 38 C.R. Cleven, 15 Y. Cobigo, 11 B.A. Cole, 10 M.P. Comets, 39 P. Constantin, 19 M. Csanad, 13 
T. Csorgo, 24 D. d'Entcrria, 10 ' 50 T. Dahms, 49 K. Das, 14 G. David, 5 H. Delagrange, 50 A. Denisov, 17 
A. Deshpande, 43,49 E.J. Desmond, 5 A. Dcvismes, 49 O. Dietzsch, 46 A. Dion, 49 J.L. Drachcnberg, 1 O. Drapier, 28 

A. Drees, 49 K.A. Drees, 5 A.K. Dubcy, 56 R. du Rictz, 32 A. Durum, 17 D. Dutta, 4 V. Dzhordzhadze, 51 
Y.V. Efrcmcnko, 38 J. Egdemir, 49 K. El Chcnawi, 54 A. Enokizono, 16 H. En'yo, 42 ' 43 B. Espagnon, 39 S. Esumi, 53 
L. Ewcll, 5 D.E. Fields, 36 ' 43 F. Fleuret, 28 S.L. Fokin, 26 B. Forestier, 31 B.D. Fox, 43 Z. Fraenkel, 56 J.E. Frantz, 10 
A. Franz, 5 A.D. Frawley, 14 Y. Fukao, 27 ' 42 S.-Y. Fung, 6 S. Gadrat, 31 S. Garpman, 32 - * F. Gastineau, 50 M. Germain, 50 

T.K. Ghosh, 54 A. Glenn, 51 G. Gogiberidze, 51 M. Gonin, 28 J. Gosset, 11 Y. Goto, 42 ' 43 R. Granier de Cassagnac, 28 
N. Grau, 19 S.V. Greene, 54 M. Grosse Perdekamp, 18 ' 43 T. Gunji, 8 W. Guryn, 5 H.-A. Gustafsson, 32 T. Hachiya, 16 ' 42 
A. HadjHenni, 50 J.S. Haggerty, 5 M.N. Hagiwara, 1 H. Hamagaki, 8 A.G. Hansen, 30 H. Harada, 16 E.P. Hartouni, 29 
K. Haruna, 16 M. Harvey, 5 E. Haslum, 32 K. Hasuko, 42 R. Hayano, 8 N. Hayashi, 42 X. He, 15 M. Heffner, 29 

T.K. Hcmmick, 49 J.M. Heuser, 42 ' 49 M. Hibino, 55 H. Hiejima, 18 J.C. Hill, 19 R. Hobbs, 36 M. Holmes, 54 
W. Holzmann, 48 K. Homma, 16 B. Hong, 25 A. Hoover, 37 T. Horaguchi, 42 ' 52 H.M. Hur, 21 T. Ichihara, 42 ' 43 
V.V. Ikonnikov, 26 K. Imai, 27 ' 42 M. Inaba, 53 D. Isenhower, 1 L. Isenhower, 1 M. Ishihara, 42 T. Isobe, 8 M. Issah, 48 
A. Isupov, 20 B.V. Jacak, 49 W.Y. Jang, 25 Y. Jeong, 22 J. Jia, 10 ' 49 J. Jin, 10 O. Jinnouchi, 42 ' 43 B.M. Johnson, 5 
S.C. Johnson, 29 K.S. Joo, 34 D. Jouan, 39 F. Kajihara, 8 ' 42 S. Kametani, 8 - 55 N. Kamihara, 42 ' 52 M. Kaneta, 43 
J.H. Kang, 57 S.S. Kapoor, 4 K. Katou, 55 T. Kawagishi, 53 A.V. Kazantsev, 26 S. Kelly, 9 ' 10 B. Khachaturov, 56 
A. Khanzadeev, 41 J. Kikuchi, 55 D.H. Kim, 34 D.J. Kim, 57 D.W. Kim, 22 E. Kim, 47 G.-B. Kim, 28 H.J. Kim, 57 
Y.-S. Kim, 21 E. Kinney, 9 W.W. Kinnison, 30 A. Kiss, 13 E. Kistenev, 5 A. Kiyomichi, 42 ' 53 K. Kiyoyama, 35 

C. Klcin-Boesing, 33 H. Kobayashi, 42 ' 43 L. Kochenda, 41 V. Kochetkov, 17 D. Koehler, 36 T. Kohama, 16 B. Komkov 41 

M. Konno, 53 M. Kopytine, 49 D. Kotchetkov, 6 A. Kozlov, 56 P.J. Kroon, 5 C.H. Kuberg, 1 ' 30 G.J. Kunde, 30 
N. Kurihara, 8 K. Kurita, 42 - 43 ' 44 Y. Kuroki, 53 M.J. Kweon, 25 Y. Kwon, 57 G.S. Kyle, 37 R. Lacey, 48 V. Ladygin, 20 
J.G. Lajoie, 19 Y. Le Bornec, 39 A. Lebedev, 19 ' 26 S. Leckey, 49 D.M. Lee, 30 M.K. Lee, 57 S. Lee, 22 M.J. Leitch, 30 
M.A.L. Leite, 46 X.H. Li, 6 H. Lim, 47 A. Litvinenko, 20 M.X. Liu, 30 Y. Liu, 39 C.F. Maguire, 54 Y.I. Makdisi, 5 
A. Malakhov, 20 M.D. Malik, 36 V.I. Manko, 26 Y. Mao, 7 ' 42 G. Martinez, 50 M.D. Marx, 49 H. Masui, 53 F. Matathias, 49 
T. Matsumoto, 8 ' 55 M.C. McCain, 18 P.L. McGaughey, 30 E. Melnikov, 17 F. Messer, 49 Y. Miake, 53 J. Milan, 48 
T.E. Miller, 54 A. Milov, 49 ' 56 S. Mioduszewski, 5 R.E. Mischke, 30 G.C. Mishra, 15 J.T. Mitchell, 5 A.K. Mohanty, 4 

D. P. Morrison, 5 J.M. Moss, 30 T.V. Moukhanova, 26 F. Miihlbacher 49 D. Mukhopadhyay, 54 ' 56 M. Muniruzzaman, 6 

J. Murata, 42 ' 43 ' 44 S. Nagamiya, 23 Y. Nagata, 53 J.L. Nagle, 9 ' 10 M. Naglis, 56 T. Nakamura, 16 B.K. Nandi, 6 
M. Nara, 53 J. Newby, 29 ' 51 M. Nguyen, 49 P. Nilsson, 32 B. Norman, 30 A.S. Nyanin, 26 J. Nystrand, 32 E. O'Brien, 5 
C.A. Ogilvie, 19 H. Ohnishi, 5 ' 42 I.D. Ojha, 3 ' 54 H. Okada, 27 ' 42 K. Okada, 42 ' 43 O.O. Omiwade, 1 M. Ono, 53 
V. Onuchin, 17 A. Oskarsson, 32 I. Otterlund, 32 K. Oyama, 8 K. Ozawa, 8 D. Pal, 54 ' 56 A.P.T. Palounek, 30 
V. Pantuev, 49 V. Papavassiliou, 37 J. Park, 47 W.J. Park, 25 A. Parmar, 36 S.F. Pate, 37 H. Pei, 19 T. Peitzmann, 33 
J.-C. Peng, 18 - 30 H. Pcreira, 11 V. Percsedov, 20 D.Yu. Pcressounko, 26 C. Pinkenburg, 5 R.P. Pisani, 5 F. Plasil, 38 
M.L. Purschke, 5 A.K. Purwar, 49 H. Qu, 15 J. Rak, 19 I. Ravinovich, 56 K.F. Read, 38 ' 51 M. Rcuter, 49 
K. Rcygers, 33 V. Riabov, 41 - 45 Y. Riabov, 41 G. Roche, 31 A. Romana, 28 M. Rosati, 19 S.S.E. Roscndahl, 32 
P. Rosnct, 31 P. Rukoyatkin, 20 V.L. Rykov, 42 S.S. Ryu, 57 M.E. Sadler, 1 B. Sahlmueller, 33 N. Saito, 27 ' 42 ' 43 
T. Sakaguchi, 8 ' 55 M. Sakai, 35 S. Sakai, 53 V. Samsonov, 41 L. Sanfratcllo, 36 R. Santo, 33 H.D. Sato, 27 - 42 S. Sato, 5 ' 23 ' 53 
S. Sawada, 23 Y. Schutz, 50 V. Semenov, 17 R. Seto, 6 D. Sharma, 56 M.R. Shaw, 1 - 30 T.K. Shea, 5 I. Shein, 17 
T.-A. Shibata, 42 ' 52 K. Shigaki, 16 - 23 T. Shiina, 30 M. Shimomura, 53 T. Shohjoh, 53 A. Sickles, 49 C.L. Silva, 46 
D. Silvcrmyr, 30 ' 32 ' 38 K.S. Sim, 25 J. Simon- Gillo, 30 CP. Singh, 3 V. Singh, 3 M. Sivertz, 5 S. Skutnik, 19 W.C. Smith, 1 
A. Soldatov, 17 R.A. Soltz, 29 W.E. Sondheim, 30 S.P. Sorensen, 51 I.V. Sourikova, 5 F. Staley, 11 P.W. Stankus, 38 



2 



E. Stenlund, 32 M. Stepanov, 37 A. Ster, 24 S.P. Stoll, 5 T. Sugitate, 16 C. Suirc, 39 J.P. Sullivan, 30 K. Syoji, 27 ' 42 
J. Sziklai, 24 T. Tabaru, 43 S. Takagi, 53 E.M. Takagui, 46 A. Taketani, 42 - 43 M. Tamai, 55 K.H. Tanaka, 23 Y. Tanaka, 35 
K. Tanida 42 ' 43 M.J. Tannenbaum, 5 A. Tarancnko, 48 P. Tarjan, 12 J.D. Tepe, 1 ' 30 T.L. Thomas, 36 M. Togawa, 27 - 42 
J. Tojo, 27 - 42 H. Torii, 27 ' 42 R.S. Towcll, 1 V-N. Tram, 28 I. Tserruya, 56 Y. Tsuchimoto, 16 - 42 H. Tsuruoka, 53 
S.K. Tuli, 3 H. Tydcsjo, 32 N. Tyurin, 17 H. Vallc, 54 H.W. van Hcckc, 30 J. Velkovska, 5 ' 49 ' 54 M. Vclkovsky, 49 
R. Vcrtesi, 12 V. Veszpremi, 12 L. Villattc, 51 A. A. Vinogradov, 26 M.A. Volkov, 26 E. Vznuzdaev 41 M. Wagner, 27 

X.R. Wang, 15 - 37 Y. Watanabc , 42 ' 43 J. Wesscls, 33 S.N. White, 5 N. Willis, 39 D. Winter, 10 F.K. Wohn, 19 
C.L. Woody, 5 M. Wysocki, 9 W. Xie, 6 - 43 Y. Yang, 7 A. Yanovich, 17 S. Yokkaichi, 42 - 43 G.R. Young, 38 I. Younus, 36 
I.E. Yushmanov, 26 W.A. Zajc, 10 - 1 O. Zaudkte, 33 C. Zhang, 10 S. Zhou, 7 S.J. Zhou, 56 J. Zimanyi, 24 and L. Zolin 20 

(PHENIX Collaboration) 

Abilene Christian University, Abilene, TX 79699, USA 
2 Institute of Physics, Academia Sinica, Taipei 11529, Taiwan 
'' Department of Physics, Banaras Hindu University, Varanasi 221005, India 
4 Bhabha Atomic Research Centre, Bombay 400 085, India 
5 Brookhaven National Laboratory, Upton, NY 11973-5000, USA 
6 University of California - Riverside, Riverside, CA 92521, USA 
China Institute of Atomic Energy (CIAE), Beijing, People's Republic of China 
8 Center for Nuclear Study, Graduate School of Science, University of Tokyo, 7-3-1 Hongo, Bunkyo, Tokyo 113-0033, Japan 

9 University of Colorado, Boulder, CO 80309 
10 Columbia University, New York, NY 10027 and Nevis Laboratories, Irvington, NY 10533, USA 
11 Dapnia, CEA Saclay, F-91191, Gif-sur-Yvette, France 
Debrecen University, H-4010 Debrecen, Egyetem ter 1, Hungary 
18 ELTE, Ebtvos Lorand University, H - 1117 Budapest, Pdzmdny P. s. 1/A, Hungary 
14 Florida State University, Tallahassee, FL 32306, USA 
15 Georgia State University, Atlanta, GA 30303, USA 
16 Hiroshima University, Kagamiyama, Higashi-Hiroshima 739-8526, Japan 
17 Institute for High Energy Physics (IHEP), Protvino, Russia 
18 University of Illinois at Urbana- Champaign, Urbana, IL 61801 
19 Iowa State University, Ames, IA 50011, USA 
20 Joint Institute for Nuclear Research, 14-1980 Dubna, Moscow Region, Russia 
21 KAERI, Cyclotron Application Laboratory, Seoul, South Korea 
22 Kangnung National University, Kangnung 210-702, South Korea 
2S KEK, High Energy Accelerator Research Organization, Tsukuba-shi, Ibaraki-ken 305-0801, Japan 
24 KFKI Research Institute for Particle and Nuclear Physics (RMKI), H-1525 Budapest 114, POBox 49, Hungary 

5 Korea University, Seoul, 136-701, Korea 
26 Russian Research Center "Kurchatov Institute", Moscow, Russia 
21 Kyoto University, Kyoto 606, Japan 
28 Laboratoire Leprince-Ringuet, Ecole Polytechnique, CNRS-IN2P3, Route de Saclay, F-91128, Palaiseau, France 
29 Lawrence Livermore National Laboratory, Livermore, CA 94550, USA 
S0 Los Alamos National Laboratory, Los Alamos, NM 87545, USA 
81 LPC, Universite Blaise Pascal, CNRS-IN2P3, Clermont-Fd, 63177 Aubiere Cedex, France 
Department of Physics, Lund University, Box 118, SE-221 00 Lund, Sweden 
38 Institut fur Kernphysik, University of Muenster, D-48149 Muenster, Germany 

84 Myongji University, Yongin, Kyonggido 449-728, Korea 
85 Nagasaki Institute of Applied Science, Nagasaki-shi, Nagasaki 851-0193, Japan 
36 University of New Mexico, Albuquerque, NM 87131, USA 
31 New Mexico State University, Las Cruces, NM 88003, USA 
88 Oak Ridge National Laboratory, Oak Ridge, TN 37831, USA 
39 IPN-Orsay, Universite Paris Sud, CNRS-IN2P3, BP1, F-91406, Orsay, France 
40 Peking University, Beijing, People's Republic of China 
4 PNPI, Petersburg Nuclear Physics Institute, Gatchina, Russia 
42 RIKEN (The Institute of Physical and Chemical Research), Wako, Saitama 351-0198, JAPAN 
48 RIKEN BNL Research Center, Brookhaven National Laboratory, Upton, NY 11973-5000, USA 
44 Physics Department, Rikkyo University, 3-34-1 Nishi-Ikebukuro, Toshima, Tokyo 171-8501, Japan 
45 St. Petersburg State Technical University, St. Petersburg, Russia 
46 Universidade de Sao Paulo, Instituto de Fisica, Caixa Postal 66318, Sao Paulo CEP05315-970, Brazil 
41 System Electronics Laboratory, Seoul National University, Seoul, South Korea 
48 Chemistry Department, Stony Brook University, Stony Brook, SUNY, NY 11794-3400, USA 
4 Department of Physics and Astronomy, Stony Brook University, SUNY, Stony Brook, NY 11794, USA 
50 SUBATECH (Ecole des Mines de Nantes, CNRS-IN2P3, Universite de Nantes) BP 20722 - 44307, Nantes, France 

51 University of Tennessee, Knoxville, TN 37996, USA 



3 



Department of Physics, Tokyo Institute of Technology, Tokyo, 152-8551, Japan 
J ' Institute of Physics, University of Tsukuba, Tsukuba, Ibaraki 305, Japan 
54 Vanderbilt University, Nashville, TN 37235, USA 
55 Waseda University, Advanced Research Institute for Science and 
Engineering, 17 Kikui-cho, Shinjuku-ku, Tokyo 162-0044, Japan 
56 Weizmann Institute, Rehovot 76100, Israel 
57 Yonsei University, IPAP, Seoul 120-749, Korea 

New measurements are presented for charged hadron azimuthal correlations at mid-rapidity in 
Au+Au collisions at ^/snn = 62.4 and 200 GeV. They are compared to earlier measurements ob- 
tained at y/SNN = 130 GeV and in Pb+Pb collisions at y^ijvjv = 17.2 GeV. Sizeable anisotropics 
are observed with centrality and transverse momentum (pr) dependence characteristic of elliptic 
flow (v2)- For a broad range of centralities, the observed magnitudes and trends of the differen- 
tial anisotropy, V2(pr), change very little over the collision energy range ^/snn = 62 - 200 GeV, 
indicating saturation of the excitation function for V2 at these energies. Such a saturation may be 
indicative of the dominance of a very soft equation of state for ^/sjvjv ~ 60 - 200 GeV. 



PACS numbers: 25.75.Ld 

Extremely high energy-density nuclear matter is pro- 
duced in energetic Au+Au collisions at the Relativistic 
Heavy Ion Collider (RHIC) [1, 2]. The dynamical evo- 
lution of this matter is predicted to reflect the presence 
and evolution of the Quark Gluon Plasma (QGP) - a 
new phase of nuclear matter [3, 4, 5]. Azimuthal correla- 
tion measurements are important in several ways. They 
serve as a "barometric sensor" for pressure gradients de- 
veloped in the collision and hence yield insight into cru- 
cial issues of thermalization and the equation of state 
(EOS) [6, 7, 8]. They provide important constraints for 
the density of the medium and the effective energy loss 
of partons which traverse it [9, 10]. They can provide 
valuable information on the gluon saturation scale in the 
nucleus [11]. 

Recent measurements at RHIC (tJsnn = 130 and 
200 GeV) indicate a mixture of (di-)jet and harmonic 
contributions to azimuthal correlations in Au+Au colli- 
sions [12, 13, 14, 15, 16]. The asymmetric (di-)jet contri- 
butions are found to be relatively small but can be sepa- 
rated; they show an increase with px and indicate strong 
suppression of away-side jet yields [15]. Significant mod- 
ifications to the away-side jet topology have also been 
reported [17]. These observations, which are particularly 
striking for very central collisions, have been interpreted 
as evidence for parton energy loss and jet quenching in 
the produced medium [3]. The harmonic contributions 
show significant strength at mid-rapidity with character- 
istic dependencies on px and centrality [12, 18, 19, 20]. 
They are typically characterized by the second order 
Fourier coefficient, V2 = ^e 42 ^ 1_ * Rp )), where <f>\ rep- 
resents the azimuthal emission angle of a charged hadron 
and 4>rp is the azimuth of the reaction plane. The brack- 
ets denote statistical averaging over particles and events. 
At low pt (pt ^2.0 GeV/c) the magnitude and trends 
of i>2 are under-predicted by hadronic cascade models 
supplemented with string dynamics [21], but are well 
reproduced by models which incorporate hydrodynamic 
flow [5, 7]. This has been interpreted as evidence for 



the production of a thermalized state of partonic mat- 
ter [3, 4, 5]. At higher px the predictions of quark co- 
alescence [22] are consistent with the data [20, 23], and 
quantitative agreement has been achieved with transport 
model calculations which incorporate large opacities [10]. 

At Super Proton Synchrotron (SPS) energies (^/sjvat ~ 
17 GeV) azimuthal correlation measurements also in- 
dicate a mixture of (di-)jet and harmonic contribu- 
tions [24, 25]. However, the observed anisotropy of the 
harmonic contribution is approximately 50% of the value 
observed at full RHIC energy (^saw = 200 GeV). There- 
fore, an important outstanding issue is the detailed be- 
havior of V2 over the range which spans SPS - RHIC ener- 
gies. In recent work, the PHOBOS collaboration has in- 
vestigated the patterns for pr-integrated vi over a broad 
range of pseudorapidities [26]. Here, we present more re- 
vealing differential measurements for Au + Au collisions 
at ytsNN = 62.4 - 200 GeV and the first excitation func- 
tion for differential i>2 which spans beam energies from 
the Alternating Gradient Synchrotron (AGS) to RHIC 
(i/^nn - 3 - 200 GeV). 

The colliding Au beams (^/snn = 62.4, 130, and 
200 GeV) used in the measurements presented here have 
been provided by RHIC in three separate experimental 
running periods (in 2000-2004) . Charged tracks were de- 
tected in the east and west central arms of PHENIX [27], 
each of which subtend 90° in azimuth <f>, and ±0.35 
units of pseudo-rapidity 77. Track reconstruction was ac- 
complished at each collision energy via pattern recogni- 
tion using a drift chamber (DC) followed by two layers 
of multi-wire proportional chambers with pad readout 
(PCI, PC3) [27]. A combinatorial Hough transform in 
the track bend plane was used for pattern recognition 
in the DC [28]. For each analysis, the collision vertex z 
along the beam direction was constrained to be within 
\z\ < 30 cm. A confirmation hit within a 2a match- 
ing window was required in PC3, located at a radius 
of 5 m, to eliminate most albedo, conversions, and de- 
cays. Particle momenta were measured with resolutions 
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Sp/p = 0.7% © 0.91% p, Sp/p = 0.6% © 3.6% p, and 
Sp/p = 0.7% © l.O%p (GcV/c) at ^^vw = 62.4, 130, 
and 200 GeV respectively, good enough to have very lit- 
tle influence, if any, on the results presented here. 

Event centralities were obtained at y/s^N = 62.4 GeV 
via a series of cuts on the analog response of the PHENIX 
beam counters (BBC). For ^/sJjn = 130 and 200 GeV, 
cuts in the space of BBC versus ZDC analog response 
were employed; they reflect percentile cuts on the total 
interaction cross section at each beam energy [29]. Esti- 
mates for the number of participant nucleons N part , were 
also made for each of these cuts following the Glauber- 
based model detailed in Ref. [29]. Systematic uncertain- 
ties associated with these determinations are estimated 
to be less than ~ 10% for central and mid-central colli- 
sions 

The differential v 2 measurements reported in this Let- 
ter have been obtained via three separate methods of 
analysis. In the first, we used the reaction plane tech- 
nique which correlates the azimuthal angles of charged 
tracks detected in the central arms with the azimuth of 
an estimated event plane $2, determined via hits in the 
North and South BBC's located at | <q |~ 3-3.9 [20]. This 
method was used for the analysis of data taken at both 
JsWn = 62.4 and 200 GeV. Corrections [20, 30] were 
applied to account for possible azimuthal distortions in 
the distribution of the estimated reaction planes. Values 
of V2 were calculated via the expression 



r 2 



(cos(2(^-$ 2 ))) 
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where the denominator represents a resolution factor 
which corrects for the difference between the estimated 
and the true azimuth of the reaction plane &np [20, 30]. 
The estimated resolution of the combined reaction plane 
from both BBC's [20] has an average of 0.33 (0.16) 
over centrality with a maximum of about 0.42 (0.19) for 
\/snn— 200 (62.4) GeV. Thus, the estimated correction 
factor, which is the inverse of the resolution for the com- 
bined reaction plane, ranges from 2.4 (5.4) to 5.0 (13). 

In the second method, a cumulant analysis was per- 
formed on data collected at ^snn = 200 and 62.4 GeV 
to obtain the anisotropy directly [31] 
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Here, the double brackets denote an average over pairs 
of particles emitted in an event followed by further aver- 
aging over events. For a detector having full azimuthal 
acceptance, the averages (e 21 ^ 1 ^ and ^e -2 ^ 2 ) vanish due 
to symmetry considerations, to give the second order cu- 
mulant estimate v 2 {2\ [31] of v 2 



v 2 {2f 
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Since PHENIX does not have full azimuthal acceptance, 



initial underestimate of the extracted anisotropy. To cor- 
rect for this underestimate, separate correction factors 
(~ 30%) were evaluated and applied for each centrality 
and pt cut, at each collision energy, following the proce- 
dures detailed in Ref. [31]. 

In the third method, we extracted the anisotropy at 
■J snn = 62.4, 130 and 200 GeV via assorted two-particle 
correlation functions [12, 20]: 
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) and (e do not vanish and this leads to an 



Here, N cor (A(j>) is the observed A<fi distribution for 
charged particle pairs selected from the same event, and 
N m i X (A(j)) is the A(f> distribution for particle pairs se- 
lected from mixed events. Mixed events were obtained 
by randomly selecting each member of a particle pair 
from different events with the same multiplicity and ver- 
tex cuts. 

To extract the anisotropy of these correlations, two cor- 
relation functions were generated for each p? and central- 
ity selection [12, 20]. For the first, charged hadron pairs 
were formed by selecting both particles from a reference 
range pr,ref, which excluded the pr range of interest (i.e. 
a reference correlation) . For the second, assorted hadron 
pairs were formed by selecting one member from the px - 
range of interest and the other from PT,ref ■ The elliptic 
flow v 2 , was obtained via the ratio A 2 ^ a / ^ A 2 ^ re f = v 2 
where A 2 ^ a and A 2 ^ re j are the anisotropics extracted from 
the assorted and reference correlation functions (respec- 
tively) with the fit function: 

C(A<j>) = ax (l + 2A 2 cos(2A0) + Ae(-°- s ( A */») a )) (4) 

where the Gaussian and harmonic terms are used to char- 
acterize the asymmetry (at small A<j>) and the anisotropy 
of the correlation function respectively [12, 14, 15]. 

Figures la - Id show representative A<j> correla- 
tion functions obtained for charged hadrons detected 
in the PHENIX central arms (-0.35 < 77 < 0.35) at 
y/sNN = 62.4 GeV. Correlation functions for mid-central 
events (centrality = 20 - 40%) are shown for hadrons 
with 0.5 < p T < 0.7 GcV/c and 1.0 < p T < 1.5 
GeV/c in Figs, la and c respectively. The same pr cuts 
have been made for the correlation functions shown in 
Figs, lb and d but for more peripheral collisions (cen- 
trality = 40 - 60%). For both sets of correlation func- 
tions 0.65 < pr,ref < 2.5 GeV/c. Figs, la - Id show 
a clear anisotropic pattern with relatively small asym- 
metries (0°/180° ratios). Such asymmetries have been 
attributed to small jet contributions to the correlation 
functions [12, 15], and are expected to decrease with de- 
creasing y/sNN- The curves in Fig. 1 indicate a fit to the 
correlation function with Eq. 4; they show an increase of 
the anisotropy with increasing impact parameter and pt ■ 
These trends are similar to those of prior AGS, SPS and 
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(a) Centrality (%) = 20-40 : (b) Centrality (%) = 40-60 
Au+Au ^1,^, = 62.4 GeV— 
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FIG. 1: Assorted-pr correlation functions (0.65 < Pr,ref < 
2.5 GeV/c) for charged hadrons of 0.5 < p T < 0.7 GeV/c 
(top panels) and 1.0 < pr < 1-5 (bottom panels) obtained in 
Au+Au collisions at ^/sjvjv = 62.4 GeV. The left and right 
panels show correlation functions for centrality cuts of 20- 
40% and 40-60% respectively. The lines represent fits to the 
correlation functions (see text). 



RHIC measurements [18, 24, 25, 32] and are consistent 
with the expected patterns for in-plane elliptic flow [5, 7]. 
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FIG. 2: Differential anisotropy V2(pr) for charged hadrons 
in Au+Au collisions at ^/sjvjv = 62.4 GeV with centrality 
cuts of 0-20% (open symbols) and 20-40% (filled symbols), 
obtained via the methods of correlation functions (CF), cu- 
mulants (Cum) and reaction plane (RP). 

Figure 2 compares the differential anisotropy V2{pt), 
obtained at y/sjyN = 62.4 GeV for all three methods of 
extraction. The error bars shown indicate statistical er- 
rors. Systematic errors are estimated to be ~ 10%, 5%, 
and 5% for extractions via the reaction plane, cumulant 
and correlation function methods of analysis respectively. 
The results, which are shown for two separate centrality 



cuts (0 - 20% and 20 - 40%) in each case, indicate an 
initial increase of z>2 with px followed by the previously 
observed plateau for p T £ 2.5 GeV/c [12, 19]. The close 
agreement of V2{pt) values obtained from the cumulant 
and correlation function methods of analysis, serve to 
confirm the reliability of these methods of extraction. On 
the other hand, the agreement between results from these 
latter methods and that obtained from the reaction plane 
method is quite striking, given the large rapidity gap (~ 3 
units) between the particles used for reaction plane deter- 
mination and the mid-rapidity particles correlated with 
this plane. It is expected that the latter correlations are 
less influenced by non-flow contributions especially for 
Pt < 2.0 GeV/c. Consequently, we attribute this agree- 
ment to the absence of strong non-flow contributions to 
the hadron correlations (for px < 2.0 GeV/c) at mid- 
rapidity. A similarly good agreement between the differ- 
ent methods of analysis was obtained for all centralities 
presented in this work. 
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FIG. 3: Differential anisotropy V2(N par t) (left) and 1^2 (pt) 
(right) for several energies as indicated. V2(pr) is shown for 
the centrality selection 13 - 26%. The CERES data are taken 
from Ref. [24]. 

Figures 3a and 3b compare the centrality and pr de- 
pendence (respectively) of the anisotropy obtained at 
several collision energies. The circles, stars and squares 
in 3a show v^i^part) for (pr) selections of 0.4, 0.75 and 
1.35 GeV/c obtained via the cumulant and correlation 
function methods of analysis. The same results obtained 
via the reaction plane method are consistent with prior 
results[20]. The open and filled symbols show measure- 
ments performed at ^/satjv = 62.4 and 130 (200) GeV as 
indicated; they show rather striking agreement between 
the magnitudes of the values obtained at all three 
collision energies. Further evidence that this agreement 
persists down to ^/snn = 62.4 GeV is given in Fig. 3b. 
Here, the open and filled circles compare the differen- 
tial anisotropy ^(pt), obtained at ^/sjvtv = 62.4 and 
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200 GcV for the 13-26% most central collisions. The 
comparison indicates little change in i>2 as the collision 
energy is raised from ^/snn = 62.4 to 200 GeV. This 
contrasts with the much lower vi values measured in 
Pb+Pb collisions (filled squares) by the CERES collabo- 
ration at y/SNN = 17.2 GeV, for the same centrality cut 
(13 - 26%) [24]. 
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FIG. 4: Differential V2 vs. y / sjvjv for charged hadrons in 
nucleus- nucleus collisions. Results are shown for the cen- 
trality cut of 13 - 26% and pr selections of 1.75 GeV/c 
(open symbols) and 0.65 GeV/c (closed symbols). The 
STAR, CERES and E895 data are taken from Refs. [19], [24] 
and [32, 33, 34] respectively. 



The tJsnn dependence of vi for charged hadrons pro- 
duced in Au+Au collisions is summarized in Fig. 4 for 
two separate pr selections (0.65 and 1.75 GeV/c) and 
centrality = 13-26%. These data are taken from the 
current measurements and earlier measurements at the 
SPS [2 [} and the AGS [32, 33, 34]. The AGS measure- 
ments [32, 33, 34] are for protons. The STAR results 
were obtained for a slightly different centrality selection 
(10-30%) [19] having essentially the same mean central- 
ity. For both px cuts, the magnitude of v-i shows a 
significant increase with collision energy (~ 50% increase 
from SPS to RHIC) up to the energy ^/sWW = 62.4 GeV. 
Thereafter, it appears to saturate for larger beam ener- 
gies. 

To summarize, we have measured differential az- 
imuthal anisotropics for charged hadrons in Au + Au 
collisions spanning the energy range ^/snn = 62.4 - 200 
GeV. Detailed comparisons of these differential measure- 
ments indicate no significant collision energy dependence 
of the anisotropy over this range. By contrast, compar- 
isons to differential measurements obtained at AGS and 
SPS energies indicate that V2 increases with collision en- 
ergy up to y/SNN = 62.4 GeV. Given the fact that the 
energy density is estimated to increase by approximately 
30% over the range -^/sTT/v = 62.4 - 200 GeV, this ap- 



parent saturation of V2 above ^/snn = 62.4 GeV may be 
indicative of the role of a rather soft equation of state. 
Such a softening could result from the production of a 
mixed phase [33] for the range y/s^N = 62.4 - 200 GeV. 

We thank the staff of the Collider-Accelerator and 
Physics Departments at BNL for their vital contri- 
butions. We acknowledge support from the Depart- 
ment of Energy and NSF (U.S.A.), MEXT and JSPS 
(Japan), CNPq and FAPESP (Brazil), NSFC (China), 
IN2P3/CNRS, and CEA (France), BMBF, DA AD, and 
AvH (Germany), OTKA (Hungary), DAE and DST (In- 
dia), ISF (Israel), KRF, CHEP, and KOSEF (Korea), 
RMIST, RAS, and RMAE (Russia), VR and KAW (Swe- 
den), U.S. CRDF for the FSU, US-Hungarian NSF- 
OTKA-MTA, and US-Israel BSF. 



* Deceased 

t PHENIX Spokesperson:zajc@nevis. columbia.edu 
[1] QM2004, J. Phys. G30, S633-S1429 (2004), (Proceedings 

to be published in J Phys. G). 
[2] J. L. Nagle (2004), nucl-ex/0405012. 
[3] M. Gyulassy et al. (2004), nucl-th/0405013. 
[4] B. Muller (2004), nucl-th/0404015. 
[5] E. V. Shuryak (2004), hep-ph/0405066. 
[6] J.-Y. Ollitrault, Phys. Rev. D46, 229 (1992). 
[7] P. F. Kolb et al., Nucl. Phys. A696, 197 (2001), hep- 

ph/0103234. 
[8] T. Hirano et al. (2004), nucl-th/0404039. 
[9] M. Gyulassy et al., Phys. Rev. Lett. 86, 2537 (2001), 

micl-th/0012092. 
[10] D. Molnar et al., Nucl. Phys. A697, 495 (2002), nucl- 

th/0104073. 

[11] Y. Kovchegov et al., Nucl. Phys. A708, 413 (2002), hep- 
ph/0203213. 

[12] N. N. Ajitanand et al., Nucl. Phys. A715, 765 (2003), 

nucl-ex/0210007. 
[13] M. Chiu et al., Nucl. Phys. A715, 761 (2003), nucl- 

ex/0211008. 

[14] C. Adler et al., Phys. Rev. Lett. 90, 032301 (2003), nucl- 
ex/0206006. 

[15] C. Adler et al., Phys. Rev. Lett. 90, 082302 (2003), nucl- 

ex/0210033. 
[16] S. S. Adler et al. (2004), nucl-ex/0408007. 
[17] J. Rak, J. Phys. G30, S1309 (2004), hep-ex/0403038. 
[18] K. Adcox et al., Phys. Rev. Lett. 89, 212301 (2002), nucl- 

ex/0204005. 

[19] C. Adler et al., Phys. Rev. C66, 034904 (2002), nucl- 
ex/0206001. 

[20] S. S. Adler et al., Phys. Rev. Lett. 91, 182301 (2003), 

nucl-ex/0305013. 
[21] M. Bleicher et al., Phys. Lett. B526, 309 (2002), hep- 

ph/0006147. 

[22] R. J. Fries et al., Phys. Rev. C68, 044902 (2003), nucl- 
th/0306027. 

[23] J. Adams et al., Phys. Rev. Lett. 92, 052302 (2004), nucl- 
ex/0306007. 

[24] G. Agakichiev et al., Phys. Rev. Lett. 92, 032301 (2004), 
nucl-ex/0303014. 



7 



[25] C. Alt et al., Phys. Rev. C68, 034903 (2003), nucl- 

ex/0303001. 
[26] B. Back et al. (2004), nucl-ex/0406021. 
[27] K. Adcox et al., Nucl. Instrum. Meth. A499, 469 (2003). 
[28] J. Mitchell et al., Nucl. Instrum. Meth. A482, 491 

(2002), nucl-ex/0201013. 
[29] K. Adcox et al, Phys. Rev. C69, 024904 (2004), nucl- 

ex/0307010. 

[30] A. M. Poskanzer et al., Phys. Rev. C58, 1671 (1998), 
nucl-ex/9805001. 



[31] N. Borghini et al., Phys. Rev. C64, 054901 (2001), nucl- 
th/0105040. 

[32] P. Chung et al., Phys. Rev. C66, 021901 (2002), nucl- 
ex/0112002. 

[33] P. Danielewicz et al., Science 298, 1592 (2002), nucl- 
th/0208016. 

[34] C. Pinkenburg et al. (E895), Phys. Rev. Lett. 83, 1295 
(1999), nucl-ex/9903010. 



